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ABSTRACT 

Ozone profiles me derived from backscattered radiances in 
the ultraviolet spectra (290-340 mn) measured by the nadir- 
viewing Global Ozone Monitoring Expe rim e n t using 
optimal estimation. Tropospheric O 3 is directly retrieved 
with die tropopause as one of the retrieval levels. To 
optimize the retrieval and improve the fitting precision 
needed for tropospheric O3, we perform extensive 
wavelength and radiometric calibrations and improve 
forward model inputs. Retrieved O 3 profiles and 
tropospheric O 3 agree well with coincident ozonesonde 
meas ur eme n t s , and the integrated total O 3 agrees very well 
with Earth Probe TOMS and Dobson/Brewer total O3. The 
global distribution of tropospheric O 3 clearly shows the 
influences of biomass burning, convection, and air 
pollution, and is generally consistent wiht our current 
understanding. 

1. INTRODUCTION 

The retrieval of O 3 profile including the troposphere from 
Global Ozone Monitoring Ex p e rim e n t (GOME) has been 
demonstrated in recent years using physics-based 
approaches [1-4]. These algorithms require very accurate 
radiometric and wavelength calibrations, and accurate 
modeling of the atmosphere other than O 3 , including clouds, 
aerosols, and temperature profiles. 

This study performs detailed wavelength and radiometric 
calibrations, improves forward model inputs using our best 
available knowledge, and derives O 3 profiles and 
tropospheric <>3 from measured GOME radiances in the 
ultraviolet We validate the retrievals against ozonesonde, 
Dobson/Brewer and TOMS m eas ur ements and present the 
global distribution of tropospheric O3. 

2. DATA AND METHODOLOGY 

Ozone profile retrieval from nadir spectra is an ill- 
conditioned problem. Therefore, we choose the well-known 
Optimal Estimation (OE) technique for inversion [5]. OE 
uses available a priori knowledge to stabilize retrievals. 
Because the information content from nadir-viewing spectra 


is limited, using proper a priori climatology is important 
We uk the TOMS V 8 climatology as a priori and a priori 
variance. This climatology is month- and latitude- 
dependent and is derived from Stratospheric Aerosol and 
Gas Experiment (SAGE), Polar Ozone and Aerosol 
Measurement (FOAM), and ozonesonde observations [ 6 ]. A 
correlation length of 6 km is used to construct the a priori 
covariance matrix. 

We improve the wavelength and radiometric calibrations as 
follows. (1) We derive variable slit widths, and shills 
betwe e n radiances/irradiaiices, at every 2 -mn region with a 
high-resolution solar reference spectrum [7]. (2) Shifts 

between trace gas absorption cross-sections and radiances 
are fitted in the retrieval. (3) We perform on-tine correction 
of die filling in of solar and telluric absorption features 
using Ring spectra calculated wife a first-order rotational 
Raman scattering model [ 8 ]. Ring spectra are updated when 
total O 3 changes by £ 20 Dobson Units (DU). (4) 

Undersampling of GOME is corrected using a high- 
resolution solar re fe rence spectrum [9]. (5) We improve the 
polarization correction to GOME measureme nts using the 
GOMECAL package [10] 

(http://www.knmi.nl/gomc fd/gomecal/ ). 

We improve ch ar ac t e riz ation of fee atmosphere wife cloud 
information from the GOME Cloud Retrieval Algorithm 
[11], monthly mean SAGE n aerosols [12] and GEOS- 
CHEM tropospheric aerosols [13]; daily European Centre 
for Medium-Range Weather Forecasts (ECMWF) 
t emperat ure profiles (http://www.ecmwf.intl for extracting 
tropospheric O 3 from fee t e m p er a t u re-dependent Huggins 
bands [14], and daily surface pressure from National Coolers 
for Environmental Prediction/National Center for 
Atmospheric Research (NCEP/NCAR) reanalysis data 
(http://www.cdc.noaa.gov) . Initial albedo is based on 
derived GOME surface albedo database [15]. Cloud 
fraction is readjusted from measured re flect a nce at 370.02 
mn where absorption is minimal. 

The retrieval uses measurements in the windows 290-307 
mn and 327-336 nm. Meas uremen t s below 290 nm are not 
used because of large measurements errors and NO emission 


lines. We also find that including band lb measurements 
between 307 and 314 nm does not improve the retrievals 
probably because of inconsistent calibration with band la 
measurements. In addition, GOME slit function and 
wavelength shifts vary rapidly at the beginning of band 2b 
(312-327 nm), which makes it difficult to use them to 
improve the retrievals. The spatial resolution of retrievals is 
960 km x 80 km. 

Partial columns are typically retrieved on an 11 -layer 
Umkehr-like grid except in the troposphere, where retrieval 
grids are modified using tropopause pressure and surface 
pressure. Besides the 0 3 variables, the state vector includes 
an albedo parameter for band lb, three albedo parameters 
for band 2b to account for its wavelength-dependence, four 
cross section shift parameters for each band, four 
wavelength shift parameters for each band, scaling and shift 
parameters for minor species (e.g., N0 2 , S0 2 , BrO), two 
parameters for undersampling correction, one Ring scaling 
parameter for band la, and three Ring scaling parameters for 
band 2b to account for multiple scattering, and three 
parameters to account for the degradation correction in band 
la. The degradation parameters are used only in band la 
because band la measurements are more severely degraded. 
The total 0 3 information from band 2b implicitly constrains 
the degradation correction in band la. 

We use LIDORT [16] to calculate radiances and weighting 
functions. The scalar radiances from LIDORT are corrected 
for neglecting polarization using a look up table. 

The retrievals are compared with correlative TOMS V8 and 
Dobson/Brewer total 0 3 measurements and ozonesonde 
observations. The equator crossing time difference between 
TOMS and GOME is approximately one hour. We use the 
gridded TOMS data (1.0° latitude x 1.25° longitude); all 
TOMS measurements (~ 20 values) within a GOME 
footprint are averaged. The coincident criteria for Dobson 
/Brewer and ozonesonde measurements are 1.5° degree in 
latitude, 600 km in longitude, and 8 hours. Ozonesonde 
measurements during 1996-2000 at 1 1 stations are used in 
these studies (Table 1). Six stations also have total 0 3 
measurements from the Dobson or Brewer measurements. 

Table 1. List of ozonesonde stations (abbreviation name), 
locations, the time period during which data are used, and the 
availability and type of total Os (TOZ) measurements. The 
measurements at Scoresbysund are obtained from Network for the 
Detection of Stratospheric Change (NDSC, 
http://www.ndsc.ncep.noaa.gov ). The measurements at Java are 
obtained from M. Fujiwa before 1998 and from Southern 
Hemisphere Additional Ozonesondes (SHADOZ) [17] since 1998. 
Other data are obtained from World Ozone and Ultraviolet 
Radiation Data Center (WOUDC, 

http://www.woudc.org/index_e.html). 


Station (abbr.) 

Location 

Time 

TOZ 

Ny Alesund (ny) 

78.9°N, 11.9*E 

96-00 

N/A 1 

Scoresbysund (sc) 

70.5°N, 22.0“W 

96-00 

N/A 

Sodankylfi(so) 

67.4°N, 26.7°E 

96-98 

Brewer 

Hohenpeifienberg (ho) 

47.9"N, 1 1.0°E 

96-00 

Dobson 

Hilo (hi) 

19.6®N,155°W 

96-00 

Dobson 2 

Nairobi (nr) 

1.3°S, 36.8°E 

97-00 

Dobson 3 

Java (ja) 

7.6*S, 112.7°E 

96-00 

N/A 

Ascension (as) 

8.0°S, 14.4*W 

97-00 

N/A 

America Samoa (sa) 

14.2°S,170.6*W 

96-00 

Dobson 

Lauder (la) 

45.0°S, 169.7*E 

96-00 

Dobson 

Neumayer (ne) 

70.T*S, 8.3°W 

96-00 

N/A 


1 . Total Oj is available only during 1996-1 997 and is not used. 

2. The total Oj is measured at a close station Mauna Loa (1 9.5*N, 


155.6*E); its surface altitude is 3.4 1cm. 

3. Total 0 3 is available only during 1 997- 1 999. 


3. RESULTS AND DISCUSSION 


3.1. Comparisons with TOMS and Dobson/Brewer Total 0 3 




Fig. 1. Comparisons of 
integrated total O 3 (purple) 
with TOMS (red) and 
Dobson total 0 3 (green) at 
Hobenpeipenberg. (a) Total 
0 3 . (b) Difference between 
GOME and TOMS, (c) 
Difference between 

GOME/TOMS and Dobson. 
The a priori 0 3 is also 
shown as yellow in (a). 


Ffg. la shows the time series of integrated total 0 3 together 
with TOMS and Dobson total 0 3 at Hobenpeipenberg. 
Their differences are shown in Fig. lb and Fig. lc. The 
time series of differences do not change with time, 
suggesting that band 2 measurements do not degrade much 
and the severe degradation in band la since 1998 [3] are 
well handled using the above-mentioned degradation 
correction scheme. We can see that our retrievals agree very 
well with both TOMS and Dobson total ozone with average 
biases less than 4 DU. Table 2 shows the biases, standard 
deviations, and correlation coefficients between our 
retrieved, TOMS, and Dobson/Brewer total 0 3 . The biases 
are within 7 DU at all stations; the correlation coefficients 
are greater than 0.92 between GOME and TOMS and 
greater than 0.87 between GOME and Dobson/Brewer. 
These biases are within the retrieval uncertainties of 
different retrievals and the spatiotemporal variability of total 
0 3 . Compared to the a priori , we can certainly see 
significant improvements, with smaller biases and much 
smaller standard deviations. The standard deviations are 



smaller for comparisons with TOMS measureme nts because 
of similar spatial domains. In the comparison with Dobson 
measurements, we are comparing area vs. point 
measurements; it is expected that the standard deviations 
will be larger. We note that die standard deviations are 
larger at higher latitudes because of larger spatiotemporal 
variability. For example in Fig. lc, there are a few large 
differe n ces, greater than 40 DU, be tw een GOME/TOMS 
retrievals and Dobson meas ure me n ts in the January and 
February of 1999. The longitude differences are g r eate r 
than 5* and the time difference are greater than 4 hours. 
There are large spatial gradients over this region for these 
days as seen from TOMS data, suggesting that these large 
differences result from the large spatiotemporal Oj 
variability at mid-latitudes in the winter and early spring. 

Table 2. Biases, standard deviations, and comfatkm co efficients 
between GOME total Oj and if prfori/TOMS/Dobsoo/Brcwer total 
0 3 . The O 3 between retrieval surface to station sirface attitude is 
taken into account in the following comparisons. The units for 
biases and standard deviations are DU. 


Statiaa 

GOME -4 
« « 

man 

GOME- 

TOMS 

GOME- 

ny 

-23.3i39J0.66 -3J±8J/0.99 

N/A 

sc 

-13.3*40*0.68 

2.0*6.8*0.99 

N/A 

so 

-29.3*41*0.53 

0.0±7.8/0.99 

-52*142/0.97 

ho 

-6.2i30.6/0.66 

2.7*53*0.99 

3.9*13.9*0.94 

hi 

-5.1*13.1*0.66 

-2.3*4.2*0.97 

0.8*63*0.95 

nr 

1 .3*122*022 

-13*4.8/0.93 

-2.9*6.0/0.87 

j» 

-1 3. 4i7. 6/0.62 

-5.6*33*0.92 

N/A/ 

as 

1 .8*8.2*0.59 

0-2*3.7*9.94 

N/A 

sa 

-9.8*8.3/0.54 

-3.9*3.4*0.93 

23*42*0.91 

la 

-8.5*23.2/0.80 

-3.6*62/0.99 

0.8*18.8/0.88 

ne 

-26.0*31/0.85 

-6.3*7.9*0.99 

N/A 


3.2. Comparisons of 0 3 profiles with ozonesonde 


Because ozonesonde measures O 3 only up to - 35 km, we 
can only compare the bottom seven layers. We in tegrate the 
ozonesonde profiles to partial columns at each layer 
corresponding to their collocated GOME retrievals. In order 
to assess the accuracy of the GOME retrievals, the 
ozonesonde profiles are not convolved with retrieval 
averaging kernels because smoothing contributes to retrieval 
errors. Fig. 2 shows the comparison of retrieved O 3 profiles 
with ozonesonde measurements at H ob e np ciPenbcig. We 
can see that the retrievals agree very well with ozonesonde 
measurements. The average biases are within 2.5 DU and 
10% at each layer (Fig. 3). 

Fig. 3 shows the average biases and standard deviations 
between retrieved/a priori profiles and ozonesonde profiles. 
For Scoresbysund, Sodankyla, Hohenpeiftenberg, and 
Lauder, the comparisons are very good with average biases 
within 6 DU and 15% at each layer. The standard 
deviations between retrievals and ozonesonde are usually 


reduced relative to those between a priori and ozonesonde 
Oy 



Fig 2. Comparisons of 
GOME retrieved (purple 
circles) and ozonesonde 
O 3 profiles (green 
triages) for the bottom 
seven layers (0—35 km) at 
Ho h ta p ci p cfl bcrg. The 0 
priori arc shown as 
yellow. 
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Fig. 3. Mean biases (solid) and standard deviations (dashed) at all 
stations for (a) GOME - sonde, (b) GOME a priori - sonde, (c) 
and (d) are similar to (a) and (b) except for relative biases and 
standard deviations. 


At Ny Alesund, there is generally good a greem ent except at 
for layer 2 (-10-15 km), where the mean bias is -13 DU, 
although die difference between a priori and ozonesonde is 
much smaller (-2.6 DU). This difference is usually larger 
during the spring than the summer, and is slightly correlated 
with solar zenith angle (R=-024) and cloud fraction (R=- 
0.34). The bias is probably due to the incorrect assignment 
of ice surface to clouds and is c ur rently under investigation. 
At Neumayer, large biases occur at layer 4 (-20-25 km) 
with a mean bias of 10 DU. The relative bias can be greater 
than 200 %, especially when there is strong O 3 depletion and 
the column O 3 is less than 10 DU. The bias is partly due to 
similar large bias in the a priori Oy For the five tropical 
stations, the agreements are within 5 DU and 10% for layers 
5 and 6 and are within 4 DU and 30% for layers 0 and 1, but 
large biases occur at layers 3, 4, and 5. Thompson et aL [17] 




also reported large biases between TOMS and SHADOZ 
integrated and evaluated total 0 3 . Currently, our integrated 
total 0 3 columns agree very well with Dobson and TOMS 
total O 3 columns. It is not clear whether these large biases 
result from retrieval errors or ozonesonde measurement 
errors. We notice that the biases for those layers at 
America Samoa and Hilo are actually much smaller before 
early 1998, when sensor solution for measurements was 
switched from 1% Kl-buffered solution to 2% KJ 
unbuffered solution. 

3.3 Comparisons of tropospheric 0 3 with ozonesonde 




Fig. 4 Comparisons between GOME retrieved (purple circles) and 
ozonesonde (green triangles) tropospheric 0 3 columns at (a) 
Hohenpeipenberg, (b) Hilo, (c) Java, (d) Ascension, and (e) 
American Samoa. The a priori are also shown as yellow. 


Fig. 4 compares the integrated tropospheric 0 3 column at 
five stations. We can see that our retrievals agree very well 
with ozonesonde measurements. Most of the small-scale 
variabilities are captured by our retrievals. For example, at 
Java, although the a priori values are usually around 30 DU, 
our retrievals successfully capture small values (~15 DU) as 
well as enhanced values (~50 DU) during the period of 
intense biomass burning resulting horn 1997-1998 El Niflo 
event. Table 3 summarizes the biases, standard deviations 
and correlation coefficients at all the stations. The average 
biases are usually within 3 DU and the standard deviations 
are within 9 DU except at Ny Alesund, Scorcsbysund, and 
Neumayer, where there are large biases. The large biases at 
those stations may result from the incorrect treatments of 
cloud and ice surface. There is a large standard deviation at 
Ascension Island, which is because the South Atlantic 
Anomaly can largely affect the retrieval. 

Table 3. Biases, standard deviations and correlation coefficients 
between retrieved/a priori and ozonesonde tropospheric 0 3 . The 
units for biases and standard deviations are DU. 


Station 

GOME -SONDE 

A Priori - SONDE 

ny 

-7.0*6.0/0.22 

0.64*6.5/0.39 

sc 

-3.5±5.5/0.46 

2.6*6.2/0.61 

so 

-2.0*5.4/0.59 

2.7*6.1/0.59 

ho 

0.0*6.8/0.57 

2.4*5.8/0.65 

hi 

-2.4*7.7/0.73 

0.5*8.0/0.68 

nr 

-3.0*8.0/0.42 

-1.8*5.4/0.29 

ja 

-1.4*4.8/0.86 

5.8*8.2/0.47 

as 

-1.3*9.3/0.52 

-7.13*7.0/0.65 

sa 

0.6*5.3/0.68 

5.0*6.3/0.45 

la 

-1.4*6.4/0.10 

-1.0*4.8/0.50 

ne 

-4.7*6.6/0.78 

2.2*4.1/0.92 


The retrieval errors in tropospheric 0 3 due to random noise 
and smoothing are less than 3 DU for tropical regions and 
less than 5 DU for mid-latitude regions. The accuracy in 
ozonesonde measurements is 5-15% [18]. The natural 
tropospheric 0 3 variability is 20-30% at mid-latitudes [19]; 
the tropospheric 0 3 can change by a factor of 3 at most 
SHADOZ stations [20]. Considering the spatial and time 
domain difference and the spatiotemporal variability of 
tropospheric 0 3 , we can say that our retrievals are very 
consistent with ozonesonde measurements within the 
measurement/retrieval uncertainties and spatiotemporal 
variability. 

3.4 Global distribution of tropospheric 0 3 

Fig. 5 (top) shows the global distribution of tropospheric 
O 3 during 2/24-26/1997. In the tropics, there is low 0 3 
over the Pacific Ocean, where there are intense convection 
activities. Over North Africa (~0°E, 10°N), a region with 
intense biomass burning during this period, there is 
enhanced tropospheric O 3 , which is not presented in most 
of the tropospheric O 3 retrievals [ 21 ]. However, the O 3 





values are smaller than those over the South Atlantic 
Ocean, consistent with the observed Atlantic tropospheric 
paradox [22]. In the Antarctic, the tropospheric O 3 is 
consistently less than 15 DU. Near 30*N and 30*S, which 
corresponds to downward motion of Hadley circulation, 
there are bands of high tropospheric O 3 even over the 
ocean. These high O 3 values are probably due to die 
complex interplay of photochemistry and 

localized/transported pollution under favorable weather 
conditions. During 9/16-18/1997, the two bends of high O 3 
and the low O 3 over Pacific Ocean with shift slightly north 
with the motion of mt ertrop ic al c o nve rgence zone (ITCZ). 
There is enhanced t ropos ph eric O 3 over Indonesia, 
consistent with the intense biomass burning and 
meteorological conditions (e.g. dry air, less p re cipi t a tion) 
caused by die 1997-1998 El Nifio event Over South 
America, South Africa and die Atlantic Ocean, high O 3 
values of 40-60 DU result from intense biomass burning 
over South Africa and South America. 



Fig. 5 Three-day composites of global distribution of tropospheric 

0 3 . (a) 2/24-26/1997. (b) 9/16-18/1997. The noisy patterns over 
the South Atlantic Ocean are caused by the South Atlantic 
Anomaly. 

4. CONCLUSION 

An algorithm was developed to retrieve O 3 profiles and 
tropospheric O 3 from GOME using the opt ima l estimation 
technique. We particularly focus on tropospheric O 3 
derivation by performing extended wavelength and 
radiometric calibrations and improving forward model 


inputs. Tropospheric 0 3 is directly retrieved by using 
tropopause as one of the retrieval layers. We validate our 
retrievals against TOMS and Dobson/Brewer total O j and 
ozonesoode measurements at 1 1 high-latitude, mid-latitude, 
and tropical stations. The integrated total ozone agrees very 
well with TOMS and Dobson/Brewer measurements with 
average biases less than 7 DU. Retrieved 03 profiles 
generally agree well with ozonesoode measurements except 
some large biases in the stratosphere at Ny Alesund, 
Neumayer, and tropical stations. The large biases at Ny 
Alesund mid Neumayer may be related to the incorrect 
assignment of snow/ice surface to clouds. It is unclear 
whether the large biases at those tropical stations are due to 
ozonesoode measurement errors or our retrieval errors. The 
retrieved tropospheric O 3 columns agree very well 
ozonesoode measurements except at Ny Alesund and 
Neumayer, cap tu rin g most of low and high Qj in the 
ozonesoode measurements. Global distributions of 
tropospheric O 3 are presented. They clearly shows signals 
due to biomass binning, convection, air pollution, and 
transport and are generally consistent with our current 
knowledge and understanding 
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